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Abstract 

We have used electronic structure calculations to investigate the 1,2-dehydration of alcohols 
as a model for water loss during the pyrolysis of carbohydrates found in biomass. Reaction 
enthalpies and energy barriers have been calculated for neat alcohols, protonated alcohols and 
alcohols complexed to alkali metal ions (Li + and Na + ). We have estimated pre-exponential A 
factors in order to obtain gas phase rate constants. For neat alcohols, the barrier to 1,2- 
dehydration is about 67 kcal mol -1 , which is consistent with the limited experimental data. 
Protonation and metal complexation significantly reduce this activation barrier and thus, 
facilitate more rapid reaction. With the addition of alkali metals, the rate of dehydration can 
increase by a factor of 10 8 while addition of a proton can lead to an increase of a factor of 10 23 . 
© 2002 Elsevier Science B.V. All rights reserved. 
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1. Introduction 

The dehydration of carbohydrates is an important process in the pyrolysis of 
biomass, though the mechanisms and kinetics of this process are not well under¬ 
stood. Loss of water from cellulose can lead to reactive functionalities that can then 
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cross-link, eventually leading to the formation of char, a more highly carbonaceous 
material that has a higher thermal stability than the original carbohydrate. 
Dehydration becomes an important consideration in biomass pyrolysis because 
char formation may be desirable or undesirable. Char pyrolysis is thought to be an 
important source of polynuclear aromatic hydrocarbons (PAHs), and releases of 
these toxins in high temperature processes are of concern. On the other hand, 
dehydration is important in the energy densificaiton of biomass. Slow roasting of 
biomass at temperature under 300 °C (torrefaction) is known to remove water from 
biomass and to increase the energy content (per gm) and the density [1,2]. This helps 
make biomass a better fuel in biomass powered plants and in co-firing with coal. 
Dehydration during torrefaction is likely due to loss of chemical water contained in 
carbohydrates. In either case, a more fundamental understanding of dehydration is 
desirable to optimize or reduce char formation. 

In spite of the importance of dehydration of carbohydrates, little mechanistic or 
kinetic information is available concerning these important reactions. The carbohy¬ 
drates in biomass contain glycol groups that exist, either on a cyclic or open 
carbohydrate. For instance, D-glucose has the following cyclic and open forms: 

CHO 



ch 2 oh 

a-pyranose form open form 


As these structures show, there are numerous adjacent alcohol groups (glycol-like 
groups). These groups can loose a water through a four-centered transition state 
involving a hydroxyl group and a hydrogen atom from neighboring carbon atom. 




+ h 2 o 


(i) 


Glycol Transition State Vinyl Alcohol and water 

There are few experimental measurements of these types of reactions, though 
similar dehydration reactions have been measured for water loss from tert -butyl 
alcohol [4] and 2,3-dimethyl-2-butanol [5]. 






-g^ch 3 

^ch 3 


,-PH 



H 3 


CH' 


+ H 2 0 


( 2 ) 


/-butanol transition state 2-methylpropene and water 

The temperature dependent rate constants, k(T ), for the dehydration of these 
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compounds have been measured and activation energies have been extracted from 
the Arrhenius equation. 


P H - 

y ? c c s^ch 3 

H 3 c' > CH 3 


H' 


.'QH 




h 3 o 

H 3 i 


i«"CH 3 

ch 3 


% 

H 3 - 


ch 3 


+ h 2 o 


(3) 


2,3-dimethylbutan-2-ol transition state 2,3-dimethylbut-2-ene and water 


k(T) = A exp 



(4) 


where A is the pre-exponential factor and E a is the activation energy. These studies 
have derived activation energies of 66.2 and 64.7 kcal mol -1 for reactions (2) and 
( 3 ). 

Some synthetic organic transformations suggest that this relatively high barrier to 
water loss may be lowered by the addition of alkali metal ions or protons. This may 
explain some of the phenomena observed with cellulose pyrolysis. For example, it 
has been demonstrated that the addition of alkali or acid to cellulose enhances the 
loss of water and the formation of char during pyrolysis [3,9,10]. It may be that these 
species lower the activation barrier for dehydration reactions such as that shown in 
(1), which lead to greater unsaturation and thus cross-linking and char formation. 
For the reaction shown in (1), the barrier to dehydration may be lowered by the 
complexation of the alkali metal ion or proton to the hydroxy group. 


+ ¥ 

h 

HO^G—Cw 

R f 


H 

r 2 


M+ 

rf H 

HO".H 

R ^R 2 


Ri 

>D=C + H 2 0*M + 
Hcf H 


( 5 ) 


However, once again, there is little experimental work to rely on to estimate the 
kinetics for these types of reactions. 

In the absence of comprehensive experimental work, we have conducted a 
molecular modeling study of the dehydration of model compounds to investigate 
reactions such as (1) and (5). We use ab initio electronic structure calculations to 
determine the structures and energetics of reactants, products and transition states 
for these important reaction pathways. In addition to molecular geometries and 
relative energies, we also obtain vibrational frequencies. These molecular parameters 
are used in statistical mechanics calculations to estimate the rate constants and 
quantify increases in the rate of dehydration upon the addition of alkali metal ions 
and protons. 
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2. Computation details 

2.1. Theoretical approach 

Electronic structure calculations were performed using the gaussian 98 [11] suite 
of programs. Two levels of theory were used in this study: (i) The hybrid density 
functional theory, B3LYP [12,13] was used initially with a triple zeta Pople type 
basis set, 6-31 lG(d, p). This technique has been shown to have an RMS error of 3.0 
kcal mol -1 for the G2 molecule set [67], and produces vibrational frequencies with a 
reported RMS error [7] of 34 cm -1 compared with a set of 122 molecules [8]. For 
molecules with few heavy atoms, B3LYP/6-31 lG(d, p) does not require significant 
computational resources and is excellent for locating minima and transition states. 
However, comparison with known experimental values indicates that B3LYP often 
predicts transition state energies lower than the experimental values [18-32] by up to 
5 kcal mol -1 : (ii) Therefore, we use a Complete Basis Set extrapolation method [14], 
CBS-QB3, to obtain more accurate transition state energies. This technique uses 
molecular geometries and frequencies obtained at the B3LYP/6-311G(d, p) level and 
extrapolates the energy to the complete basis set limit for a CCSD(T) calculation. As 
a result, this method requires considerably more computational time and memory. 
CBS-QB3 has a reported RMS error of 2.1 kcal mol -1 as determined by comparison 
to molecules in the G2 set [14]. We will show that this technique is equally accurate 
for predicting transition state energies. 

Since the size of some of the molecules in this study, significant computational 
resources were required, even for the B3LYP calculations. We used an HP Convex 
and an SGI super cluster, which were available at the National Center for 
Supercomputer Applications (NCSA). 

2.2. Reaction energies 

In our study, we locate molecular geometries of the reactants, transition states and 
products by finding stationary points on the calculated potential energy surface. A 
subsequent calculation of vibrational frequencies, obtained from second derivatives 
of the potential energy, helps ascertain the validity of the optimized geometries. 
Reactants and products will have all positive vibrational frequencies while transition 
states will have exactly one negative (imaginary) vibrational frequency, which 
corresponds to the reaction coordinate. The validity of the transition states is 
checked with two techniques. We visually inspect the atomic motions of the negative 
frequency to ensure that they correspond to the reaction coordinate and we conduct 
Intrinsic Reaction Coordinate [15,16] (IRC) calculations to determine if the 
transition state connects to the desired products and reactants on the potential 
energy surface. 

Whenever possible, we test the validity of our computational results and the 
accuracy of the techniques used by comparison to experimental values. As an 
example, consider the dehydration of r-butanol and 2,3-dimethylbutan-2-ol, 
reactions (2) and (3). We optimized the structure of the reactants and products 
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and obtained energies and enthalpies for these species. From our calculations we 
obtain the relative enthalpies between the reactants and the products or the enthalpy 
of reaction A react 7/ 2 98- Our calculated reaction enthalpies for reaction (2) and (3) are 
shown in Table 1, which also lists the experimental values for these reactions. As can 
be seen both B3LYP/6-311G(d, p) and CBS-QB3 produce reaction enthalpies that 
agree to within 3 kcal mol -1 of established thermodynamic values. 

2.3. Activation energies 

The search for transition states for reactions (2) and (3) was started using the 
geometry of transition state found for the dehydration of ethanol (see below). The 
arrangement of the atoms involved in the dehydration was similar in the transition 
states for reactions (2) and (3) as those in the dehydration of ethanol. Starting with 
the geometry of the transition state for ethanol saved a significant amount of 
computational time. We located the transition states using either Berny optimization 
[33] or an eigenvector-following algorithm [34-36]. A comparison of the important 
bond angles and lengths for the transition states obtained in this study are shown in 
Table 2. This table also shows the dihedral angle that is related to the out-of-plane 
angle of the resulting alkene structure. This will be discussed in more detail below. 

In order to estimate reaction kinetics, transition state energies are needed, and 
Table 1 also lists the calculated relative energies at 0 K, AEq, for the transition states 
of these reactions. The AE 0 value contains the zero-point energy as determined from 
the calculated vibrational frequencies using a scaling factor [14] of 0.99. The 
Arrhenius activation energy, E a in Eq. (4), is closely approximated by A E 0 . As a 
test of our computational techniques, we compared our calculated A E 0 with 
the experimental E a values. These are also shown in Table 1 and as can be 
seen, the CBS-QB3 technique predicts these activation energies within 2 kcal mol -1 , 
but the B3LYP technique underestimates these energies by about 5 kcal mol -1 . 
This is consistent with other calculations of activation energies using B3LYP 
[18-32] and the accuracy of CBS-QB3 technique is consistent with our earlier study 
[32]. 

2.4. Estimates of rate constants 

Rate constants for the dehydration of alcohols can be estimated using Transition 
State Theory (TST). The activation energy for those reactions that have an apparent 
barrier is approximated by A E 0 and the Arrhenius pre-exponential factor in Eq. (4) is 
approximated by [62] 



( 6 ) 


where o is the degeneracy (e.g. a = 3 for ethanol because of the three equivalent 
hydrogen), k B is Boltzmann’s constant, h is Plank’s constant, R is the gas constant, 
An is the change in number of molecules between the transition state and the 
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Table 1 

Calculated reaction enthalpies, A react // 2 98 5 and relative energies, AE 0 a , for transition states 


Reaction 


B3LYP b 

CBS-QB3 

Exptl 

References 



(kcal mol - 

l ) (kcal mol “ 

! ) (kcal mol ~~*) 


Neat alcohols 
(CH 3 ) 3 COH^ 






(CH 3 ) 2 CCH 2 +H 2 0 (A react tf 298 ) 

( 2 ) 

13.6 

13.7 

12.6+0.3 

[35] 

TS2 (AE 0 ) 

(CH 3 ) 2 COHCH(CH 3 ) 2 -> 


59.6 

65.9 

66.2 

[4] 

(CH 3 ) 2 CC(CH 3 ) 2 +H 2 0 (A react tf 298 ) 

( 3 ) 

10.8 

13.0 

10.7+0.4 

[39] 

TS3 (AE 0 ) 

CH 3 CH 2 OH^ 


60.4 

66.8 

64.7 

[5] 

CH 2 CH 2 +H 2 0 (A react // 298 ) 

( 8 ) 

13.5 

11.6 

11.0 

[40] 

TS8 (AE 0 ) 

(CH 3 ) 2 CHOH-+ 


62.7 

67.4 



CH 2 CHCH 3 +H 2 0 (A react // 298 ) 

( 16 ) 

14.0 

13.3 

12.2 

[43] 

TS16 (A£ 0 ) 

HOCH 2 CH 2 OH^ 


61.5 

67.0 



CH 2 CH0H+H 2 0 (A react TT 298 ) 

( 21 ) 

8.5 

7.9 

6.7 ±2.1 

[46] 

TS21 (AE 0 ) 

Levoglucosan -> 


63.4 

69.6 



1+H 2 0 (A react /f 298 ) 

( 25 ) 

13.4 




TS27 (AE 0 ) 


68.3 




Alkali metal complexes 

CH 3 CH 2 OH Li + 






CH 2 CH 2 • (H 2 0 ■ Li + ) (A react // 298 ) 


8.9 

9.9 



CH 2 CH 2 + H 2 0 • Li + (A react i/ 298 ) 

( 10 ) 

17.1 

17.5 



TS10 (A E 0 ) 


48.2 

55.3 



CH 3 CH 2 + +LiOH (A react i/ 298 ) 
CH 3 CH 2 OH ■ Na + ->• 

( 12 ) 

92.8 

94.5 



CH 2 CH 2 • (H 2 0 ■ Na + ) (A react tf 298 ) 


8.0 

9.7 



CH 2 CH 2 +H 2 0 • Na + (A react tf 298 ) 

( 11 ) 

15.2 

18.1 



TS11 (A£ 0 ) 


53.1 

60.0 



CH 3 CH 2 + +NaOH (A react // 298 ) 

(CH 3 ) 2 CHOH • Li + -> 

( 13 ) 

110.3 

114.1 



CH 3 CHCH 2 (H 2 0 ■ Li + ) (A react tf 298 ) 


9.3 

11.6 



CH 3 CHCH 2 +H 2 0 • Li + (A react tf 298 ) 

( 17 ) 

19.0 

20.8 



TS17 (A E 0 ) 


40.4 

49.5 



CH 3 CHCH 3 + +LiOH (A react i/ 298 ) 

HOCH 2 CH 2 OH Li + ->• 

( 18 ) 

75.3 

82.7 



CH 2 CHOH ■ (H 2 0 • Li + ) (A react tf 298 ) 


1.3 

3.4 



CH 2 CH0H+H 2 0 Li + (A react tf 298 ) 

( 22 ) 

31.9 

30.7 



TS22 (A£ 0 ) 


69.1 

74.2 



CH 3 CHOH + + LiOH (A react tf 298 ) 

( 23 ) 

72.1 

74.5 



Protonated alcohols 

CH 3 CH 2 OH 2 + 






CH 2 CH 2 H 3 0 + (A react tf 298 ) 


10.8 

12.1 



CH 2 CH 2 +H 3 0 + (A react // 298 ) 

( 14 ) 

32.5 

32.1 

31.6 

[49] 

TS14(A£ 0 ) 


18.6 

20.7 



CH 3 CH 2 + +H 2 0 (A react tf 298 ) 

( 15 ) 

38.5 

34.2 

33.7 

[50] 
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Table 1 ( Continued) 


Reaction 

B3LYP b 

CBS-QB3 

Exptl 

References 


(kcal mol “ 

l ) (kcal mol “ 

l ) (kcal mol ~~ 1 

) 

(CH 3 ) 2 CHOH 2 + - 





CH 3 CHCH 2 H 3 0 + (A rcact tf 298 ) 

10.3 

14.0 



CH 3 CHCH 2 +H 3 0 + (A react i/ 298 ) 

( 19 ) 37.6 

38.1 

36.7 

[53] 

TS19 (?£■()) 

7.7 

11.8 



CH 3 CHCH 3 + H 2 0 (A km H 29S ) 

7.7 

13.5 



CH 3 CHCH 3 + +H 2 0 (A rcact tf 298 ) 

HOCH 2 CH 2 OH h + -► 

( 20 ) 24.3 

25.1 

22.9 

[54] 

CH 3 CHOH + H 2 0 (A react tf 298 ) 

-13.7 

-11.4 



CH 3 CHOH + +H 2 0 (A react tf 298 ) 

( 24 ) 3.5 

1.6 

7.0 

[56] 

TS24 

24.4 

27.6 




a The relative energies for the transition states contain the zero point energy. 
b Using a 6-311G(d, p) basis set. 


reactants (in this case An = 0) and AS* is the difference in entropy between the 
transition state and the reactant. We use the entropy obtained from electronic 
structure calculations using gaussian 98. The classical approximation is used to 
obtain entropies form vibrational frequencies and molecular structures using the 
rigid rotor and harmonic oscillator (HO) approximation. The vibrational frequencies 
are used without scaling. We have found that scaling the vibrational frequencies 
changes AS* by much less than 1 cal mol -1 K -1 and has little effect upon the A 
factor. 

To properly treat aliphatic alcohols, one must consider the hindered internal 
rotors when calculating entropies. Calculations with electronic structure techniques 
usually obtain vibrational frequencies by determining second order derivatives of the 
potential at a potential minimum. Therefore, the actual potentials for the hindered 
rotors are greatly oversimplified. These potentials should have relatively low barriers 
(1-5 kcal mol -l ) as a function of rotations about C-C and C-O bonds. Vibrational 
frequencies obtained from the second derivatives at potential minima will be 
harmonic frequencies. Using the HO approximation will result in an underestimate 
of the density of states and the entropy for these molecules. 

Since we are only interested in the relative entropy, AS* in Eq. (6), for calculating 
rate constants, we can make some simplifications when treating the hindered rotors. 
For these 1,2-dehydration reactions, there are two rotors that will be greatly changed 
between the alcohol and the transition state; the rotation about the C-C bond 
involving the OH group and the rotation about the C-O bond. For ethanol, these 
are shown below. 





o 


Table 2 


(“1 

.<5—H 


R l"'.cl_ cL"'" R 3 

R 2 c^ P_ ^R 4 

Transition state geometries 


Alcohol 

M + 

r (O-Hp) (A) 

r (Cp-H p ) (A) 

r (Cp-Q) (A) 

r (Q-O) (A) 

a(C a -0-Hp) O 

Dihedral R 3 -C a -C p -R 4 (°) 

(CH 3 ) 3 COH 


1.32 

1.34 

1.43 

2.02 

59.3 

159.8 

(CH 3 ) 2 COHCH(CH 3 ) 2 


1.35 

1.31 

1.44 

2.08 

57.2 

162.9 

ch 3 ch 2 oh 


1.26 

1.43 

1.42 

1.86 

63.8 

140.3 


H + 

1.97 

1.18 

1.40 

2.53 

41.5 

175.8 


Na + 

1.54 

1.25 

1.40 

2.24 

48.8 

176.1 


Li + 

1.63 

1.22 

1.40 

2.32 

46.0 

178.9 

(CH 3 ) 2 CHOH 


1.29 

1.38 

1.42 

1.94 

61.3 

159.4 


H + 

1.99 

1.14 

1.42 

2.84 

41.1 

177.9 


Li + 

1.71 

1.18 

1.41 

2.51 

45.6 

178.0 

hoch 2 ch 2 oh 


1.25 

1.45 

1.42 

1.87 

63.1 

156.4 


H + 

2.18 

1.17 

1.41 

2.34 

44.6 

178.9 


Li + 

1.80 

1.22 

1.40 

2.27 

43.7 

175.0 


From B3LYP/6-311G(d, p) calculations. 
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A ^ 

( , C<™»h (T) 

V/ H 

Hindered rotors for ethanol 

The HO approximation for these two modes will underestimate the entropy by 
about 1 cal mol -1 K -1 for each internal rotor [68]. This will lower the pre¬ 
exponential factor by about a factor of 3. We have used the technique of Ayala et al. 
[38] that is implemented in gaussian 98 to estimate the contribution of the rotors 
similar to those in (7) to A S* and the Arrhenius pre-exponential factor, A. For the 
molecules considered here, we assume that all other internal rotations are cancelled 
out when calculating AS*. The A factors calculated with this hindered rotor 
approximation are shown in Table 3. The A factors obtained for the dehydration of 
neat 2,3-dimethylbutanol and t -butanol (6.6 x 10 13 and 1.3 x 10 14 s -1 ) show good 
agreement with the experimental values (4.57 x 10 13 and 3.98 x 10 14 s _1 ). 


3. Results and discussion 


3.1. Ethanol 


Ethanol is the simplest alcohol capable of undergoing a 1,2-dehydration reaction. 
It will, therefore, serve as a good model system for our investigation. This reaction 
proceeds through a four-centered transition state and produces ethylene and water. 


^ P H . 

Hy c G» h 

K H 


,'QH 

/ 

H*'' ''^H 




ch 2 ch 2 + H 2 o 


( 8 ) 


Ethanol TS8 Ethylene and Water 

As discussed above, we obtained stationary points for ethanol, ethylene and water. 
The calculated enthalpies for these species are used to estimate the ethalpy for 
reaction (8), A TQact H 2 9 s. The values obtained using both B3LYP and CBS-QB3 are 
listed in Table 1. The CBS-QB3 value is in agreement to better than 1 kcal mol -1 
with established thermochemistry. 

We obtain a transition state (TS8) for reaction (8) using techniques described 
above and our structure is similar to that obtained earlier using MP2 calculations [6]. 
Frequency calculations at the transition state produced exactly one negative 
vibrational frequency, 1948i cm -1 . Fig. 1 shows a scaled representation of the 
transition state and the vectors representing the atomic motion. These vectors are 
consistent with a concerted a 1,2-dehydration reaction. The atomic motion of the 
reaction coordinate is similar to that found with MP2 for 1,2-dehdyration in ethanol 
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Table 3 

Rate constants calculated using transition state theory 


Reaction 


E a (kcal 
mor’f 

A (s' 1 ) 

fc(250 °C) 
(s-‘) 

fc(500 °C) 

(s l ) 

(CH 3 ) 2 COHCH(CH 3 ) 2 -> 

(2) 

66.8 

6.6 x 10 13 



(CH 3 ) 2 CC(CH 3 ) 2 +H 2 0 






(CH 3 ) 3 COH (CH 3 ) 2 CCH 2 +H 2 0 

(3) 

65.9 

1.3 x 10 14 



ch 3 ch 2 oh -► ch 2 ch 2 +h 2 o 

(8) 

67.4 

2.7 x 10 13 

1.8 x 10- 15 

2.4 x 10 -6 

CH 3 CH 2 OHLi + ^CH 2 CH 2 + 

(10) 

55.3 

4.5 x 10 13 

3.5 x 10~ 10 

1.0 x 10 -2 

H 2 0 ■ Li + 






CH 3 CH 2 OH • Na + ^CH 2 CH 2 + 

(11) 

60.0 

6.8 x 10 12 

5.6 x 10“ 13 

7.2 x 10“ 5 

H 2 ONa + 






ch 3 ch 2 oh 2 + ^ch 2 ch 2 +h 3 o + 

(14) 

30.7 

2.8 x 10 14 

4.1 x 10 1 

5.8 x 10 5 

CH 3 CH 2 OH Li + -»CH 3 CH 2 + + 

(12) 

93.5 

10 16 

8.2 x 10 -24 

3.6 x 10~ n 

LiOH 






CH 3 CH 2 OH • Na + ->CH 3 CH 2 + + 

(13) 

112.9 

10 16 

6.3 x 10“ 32 

3.9 x 10“ 16 

NaOH 






ch 3 ch 2 oh 2 + ^ch 3 ch 2 + +h 2 o 

(15) 

32.8 

10 16 

2.0 x 10 2 

5.4 x 10 6 

(CH 3 ) 2 CHOH -> ch 3 chch 2 +h 2 o 

(16) 

67.0 

8.2 x 10 13 

8.4 x 10- 15 

9.5 x 10^ 6 

(CH 3 ) 2 CHOH • Li + ->CH 3 CHCH 2 + 

(17) 

49.5 

2.2 x 10 14 

4.6 x 10 -7 

2.3 x 10° 

H 2 OLi + 






(CH 3 ) 2 CHOH 2 + -»ch 3 chch 2 + 

(19) 

36.9 

7.3 x 10 14 

2.8 x 10 _1 

2.7 x 10 4 

h 3 o + 






(CH 3 ) 2 CHOH • Li + ->CH 3 CHCH 3 + + 

(18) 

81.5 

10 16 

9.1 x 10“ 19 

3.4 x 10“ 8 

LiOH 






(CH 3 ) 2 CHOH 2 + ^ch 3 chch 3 + +h 2 o 

(20) 

22.7 

10 16 

3.4 x 10 6 

3.9 x 10 9 

hoch 2 ch 2 oh -► ch 2 choh+h 2 o 

(21) 

69.6 

1.2 x 10 14 

1.0 x 10- 15 

2.6 x 10^ 6 

HOCH 2 CH 2 OH Li + ^CH 2 CHOH + 

(22) 

74.2 

8.5 x 10 14 

8.4 x 10“ 17 

8.9 x 10“ 7 

H 2 0 ■ Li + 






HOCH 2 CH 2 OH H + ->• 

(24) 

27.6 

9.6 x 10 14 

2.7 x 10 3 

1.5 x 10 7 

ch 3 choh + +h 2 o 






HOCH 2 CH 2 OH Li + -► 

(23) 

73.2 

10 16 

2.6 x 10“ 15 

2.0 x 10“ 5 

CH 3 CHOH + + LiOH 







a Derived from CBS-QB3 calculations. See text. 


[6]. This transition state was shown to connect ethanol with ethylene and water using 
IRC calculations. 

The calculated energy of this transition state gives a value for AE 0 of 67.4 kcal 
mol -l . Our computed barrier is consistent with: (i) the MP2 computed values [6] (67 
kcal mol -1 ) and (ii) the experimental barriers for 1,2-dehydration of £-butanol and 
2,3-butan-2-ol. The high energy of this transition state is also consistent with its 
calculated structure. Critical bond lengths and angles for the transition state in 
reaction (8) are collected in Table 2. The long C-O bond length (1.86 A) suggests 
that to a significant extent this bond has been broken, thus requiring a significant 
portion of the C-O bond dissociation enthalpy for ethanol, D// 298 (CH 3 CH 2 - 
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Fig. 1. A scaled drawing of the transition state for 1,2-dehydration of ethanol, TS8, determined by 
B3LYP/6-311G(d, p) calculations. The arrows show the atomic motions of the reaction coordinate. The 
long arrow from the reactive H atom has been reduced by 50%. 


OH) =93.9+ 0.4 kcal mol -1 . Some of this energy is offset due to partial bond 
formation between the H and O atoms. 

The other critical bond lengths and angles of the transition state presented in 
Table 2 provide additional information about the development of the 1,2- 
dehydration. The HCCH dihedral angle (defined in the table) is about 140°, which 
is between the dihedral angle in ethanol (118°) and ethylene (180°), and the C-C 
bond length (1.42 A) is between the C-C length for ethanol (1.52 A) and ethylene 
(1.33 A). These values show the progress in reaction (8) in converting ethanol to 
ethylene plus water. 

In order to estimate the rate constant for 1,2-dehydration of ethanol, we use (6) to 
estimate the pre-exponential A factor. We used the hindered rotor calculation 
available in gaussian 98 to account for the two internal rotations shown in (7). This 
amounted to a correction of 0.8 cal mol -1 K -1 in Aand a correction of less than a 
factor of 2 in the pre-exponential A factor. The calculated A factor for 1,2- 
dehydration of ethanol (Fig. 3) is 2.7 x 10 13 s -1 and with the calculate activation 
energy, AF , 0 ( — E a ) = 67.4 kcal mol+ 1 we can estimate the rate constant for this 
process. In Table 3, we present the rate constant at temperatures typical for 
torrefaction, k( 250 °C) = 1.8 x 10 -15 s -1 , and pyrolysis, £(500 °C) = 2.4 x 10 -6 
s -1 . These rate constants are small, as would be expected for reaction with such a 
high activation energy, and this suggests that this type of dehydration is inefficient. 


3.2. Stabilization of TS8 

The Mulliken charges calculated from the CCSD(T) portion of the CBS 
calculation for ethanol and the dehydration transition state (TS8) are shown below. 
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+0.46 

*H 

+0.21 r/ 

/ -0.64 

+0.21^ -0.64 / 

/ /C ~ C ^ H+0 ' 18 
H +0.20 H+o.18 


Mulliken charges for ethanol 



TS8 


( 9 ) 


As can be seen there is appreciable build up of negative charge on the oxygen atom 
in the transition state relative to ethanol. Thus, one would expect that the transition 
state would be stabilized relative to ethanol by complexation with a positive ion. We 
discuss several examples below. 


3.3. Ethanol with alkali metal ions 


As mentioned above, addition of alkali metal ions to cellulose has a profound 
effect upon the pyrolysis of cellulose. We have investigated this by considering the 
effect of added alkali cations on the reaction energies and barrier heights for 1,2- 
dehydration. The influence of Li + and Na + on the 1,2-dehydration of ethanol is 
shown in Table 1. With both cations the final products are ethylene and the cation 
complexed to a water molecule. Along the reaction coordinate to forming these 
products, a complex is formed between the hydrated alkali cation and ethylene. 
These reactions paths are shown in (10) and (11). 



.A>H 


Na + 

pH n * 


~ H 


v i+ t 


' —- CH 2 CH 2 + H 2 0*Li + 




'? H 

a 


• ^>e=i=C<y—- CH 2 CH 2 + H 2 ONa + 


alkali TS hydrated alkali 

complex complex 


( 10 ) 


(ii) 


The calculated complexation enthalpies for these alkali cations with ethanol and 
water are shown in Table 4. Since the enthalpy of complexation of the alkali with 
water is about 6 kcal mol -1 smaller than that for ethanol, the enthalpy of reactions 
(10) and (11) are more endothermic than the dehydration of neat ethanol, reaction 
(8). The effects of the alkali metal upon the structure of ethanol and water in the 
complexes is shown in Fig. 2. As can be seen the effect is slight for these alkali 
cations. This Figure also shows that the bond between the cation and the oxygen 
atom is significantly shorter (0.38 A) for Li + than for Na + , which explains the 
higher binding energy for Li + . 
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Table 4 

Calculated binding enthalpies, A bind // 298 , for complexes 


Molecule/M + 

H + a 

Li + 

Na + 

CH 3 CH 2 OH 

182.9, experimental: 185.6 b 

38.6 

26.7 

(CH 3 )CHOH 

187.3, experimental: 189.5 b 

40.3 


hoch 2 ch 2 oh 

187.1, experimental: 195.0 b 

58.4 


h 2 o 

162.6, experimental: 165.0 b 

33.1 

20.9 


Using CBS-QB3. The binding enthalpy, A bind H 2 9 s, is defined as the enthalpy for the reaction ROH- 
M + xROH + M + . 
a Defined as proton affinities. 
b [ 17 ]. 


1.81 (Li + ) 
107.5 2.19 (Na + ) 

108.1 (Li + ) 0.975 (1C) 
107.5 (Na + ) 

107.8 (H+) 




h,g- 


1.52 

1.51 (Li + ) 
1.51(Na + ) 
1.50 (H + ) 


n- 



1.43 

1.43 (Li + ) 
1.46 (Na + ) 
1.57 (H + ) 


0.98 

0.96 (Li + ) 
0.96 (Na + ) 
0.96 (H + ) 

-H 


108.1 

109.9 (Li + ) 
108.8 (Na + ) 
115.0 (If) 


1.83 (Li + ) 
2.20 (Na + ) 
0.98 (H + ) 


M 


0.97 

0.96 (Li + ) 
0.97 (Na + ) 
0.98 (H + ) 


H 105.8 H 
103.8 (Li + ) 
105.1 (Na + ) 
113.4 (H + ) 


Fig. 2. Critical bond lengths (in A) and angles (in °) for ethanol and water neat and complexed with 
M + =Li + , Na + , and H + , calculated using B3LYP/6-311G(d, p). 


In addition to the reactions shown in (10) and (11), the ethanol/alkali complex can 
also form ethyl cation [63] and hydroxylated alkali metal. 


Li + 



(12) 

—- CH 3 CH 2 + + LiOH 

H N> 

^ h 

H j^e— C™ —- CH 3 CH 2 + + NaOH 
\f H 

(13) 

The calculated enthalpies of these reaction are high (A react // 2 98 

= 94.5 kcal mol -1 


for Li + and 114.1 lccal mol -1 for Na + ) as is shown in Table 1. We could not locate 
complexes of the carbocation with the hydroxylated alkali metal. 

The effect of these alkali metal cations upon the transition state for 1,2- 
dehydration is profound. The decrease in the barrier to dehydration is shown 
graphically in Fig. 3 for ethanol complexed to Li + . For the complex with Li + , the 
barrier for dehydration is about 12 kcal mol -1 lower than for neat ethanol, while for 
Na + the barrier is lowered by about 7 kcal mol -1 . As discussed above, this trend is 
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Fig. 3. Schematic diagram of the energetics of the 1,2-dehydration of ethanol. The solid line is for neat 
ethanol, reaction (8), the dashed-dotted line is for lithiated ethanol, reaction (10), and the dashed line is 
for protonated ethanol, reactions (14) and (15). The calculated relative energies, AE 0 , are obtained using 
CBS-QB3 and include zero point energy. 


likely due to the long Na + -O bond length (2.14 A) in the transition state relative to 
the short Li + -0 bond length (1.75 A). These alkali metal cations also have the 
expected effect of making the transition state appear more like the products as is 
shown in Table 2. Note that the dihedral angles (e.g. 179° for Li + ) are now close to 
the value for ethylene (180°), though to a lesser extent for Na + than Li + . As 
described above, we have used our calculated activation energies and entropies to 
calculate rate constants for 1,2-dehydration of the alkali complexes and the results 
are shown in Table 3. The rate constant increases by roughly 5 orders of magnitude 
with Li + and 2 orders of magnitude for Na + . For K + , an important alkali metal in 
biomass, one can anticipate that the effect will be further diminished, because of its 
larger size. These increases in the rate constants show that 1,2-dehydration is 
significantly more efficient with alkali metal cations present. 

We also estimated rate constants for reactions (12) and (13) even though we could 
not locate a transition state. Since this is a simple bond dissociation, we used the 
energy of these reactions (93.5 and 112.9 kcal mol -1 ) for E a and assumed a high pre¬ 
exponential A factor (10 16 s -1 ), as one would expect for a loose transition state. 
Even with these large A factors, the rate constants for these reactions are extremely 
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small. They are 9-19 orders of magnitude smaller than the rate constant for 
reactions (10) and (11), and thus, these channels are likely to be unimportant. 

3.4. Protonated ethanol 

Protonation of ethanol produces a larger effect upon the dehydration reaction. 
The 1,2-dehydration reaction in the presence of a proton is shown in reaction (14). 



Protonated Transition Cluster Ethylene and 

Ethanol State Protonated Water 

The proton affinities (defined [17] as the enthalpy of adding a proton at 298 K) for 
ethanol and water are presented in Table 4. Our calculated proton affinities are 
compared with literature values in this table and for ethanol and water our 
calculated values are within 3 kcal mol -1 of the experimental values. This 
discrepancy between the CBS-QB3 results and experimental measurements may be 
due to the lack of diffuse functions in the basis set used [37]. The difference in proton 
affinity between water and ethanol is high, 20 kcal mol -1 , and as a result reaction 
(14) is 20 kcal mol -1 more endothermic than dehydration of neat ethanol (8). Our 
calculated reaction energy for (14), 30.7 kcal mol -1 , is consistent with a value 
obtained using [64] MP4/6-31 lG(d, p)//MP2/6-31G(d), 31.1 kcal mol -1 , and our 
energy for the cluster (11.4 kcal mol -1 ) is also in agreement with these MP4 
calculations (12.7 kcal mol -1 ). 

As can be seen in Fig. 2, the effect of an added proton upon the structures of 
ethanol and water is also profound. In ethanol the C-O bond is lengthened from 
1.43 A for neat ethanol to 1.57 A for protonated ethanol. MP4 calculations [64] 
yielded a similar value, 1.55 A. This is presumably due to the strong bond formed 
between the proton and the oxygen atom. 

As shown in Fig. 3 and Table 1, protonation of ethanol produces a larger 
reduction in barrier to dehydration relative to the addition of Li + . The barrier is 
nearly 47 kcal mol -1 lower than for neat ethanol and roughly 30 kcal mol -1 lower 
than dehydration of ethanol facilitated by Li + . Our CBS energy for this transition 
state, 20.7 kcal mol -1 , is comparable to the value obtained with MP4 calculations 
[64]. The greater energetic stability of the protonated transition state is also reflected 
in its structure. As described for neat ethanol, the C-O bond is lengthened 
considerably in the transition state (1.86 A for the TS compared with 1.43 A for 
ethanol) and because the charge separation of the transition state is stabilized by the 
proton, this effect is enhanced in the protonated transition state. The structure we 
obtain for the transition state is completely consistent with that obtained using MP4 
[64]. The C-O bond length in the TS is now 2.53 A (2.46 A using MP4 [64]). The 
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critical dihedral angle in the transition state (177.8°, 117.2 ° using MP4 [64]) is also 
significantly larger than for neat ethanol (140.3°). This value is closer to that of 
ethylene (180°) than ethanol (118°), suggesting a late transition state, and consistent 
with a highly endothermic process. 

A consequence of the stabilization of the transition state in reaction (14) is that its 
energy (20.7 kcal mol -1 ) is now lower than the energy of the isolated products (30.7 
kcal mol -1 ). As a result, calculating the rate constant for reaction (14) using the 
energy of the transition state as the activation energy (E a ) will greatly overestimate 
the rate constant for this reaction. Proper treatment of this reaction requires RRKM 
analysis. However, for a rough estimate of the rate constant we will simply use the 
reaction energy as E a in Eq. (4), and we use the TST, Eq. (6), to calculate the A 
factor. That is, A S* will still be calculated using the structure and frequencies of the 
transition state. The rate constants for this reaction (Table 3) are 11-16 orders of 
magnitude greater than neat ethanol, making dehydration of protonated ethanol 
much more efficient. Despite the approximations made in determining these rate 
constants, the increase in dehydration rate upon protonation of ethanol is quite 
dramatic. The lifetime of these reactions (the inverse of the rate constants, t 2 50 °c = 
2 x 10 - 2 s, t 5 oo °c = 2 x 10 - 6 s) is such that these reactions would be complete 
within the typical residence time of an industrial reactor (e.g. 0.1-1 s). 

In addition to loss of H 3 0 + , protonated ethanol can undergo a loss of H 2 0 to 
form a carbocation. 

H OH 2 + 

h,J^cL, h —►ch 3 ch 2 + +h 2 o (15) 

Yl 

We conducted potential energy surface calculations, where the C-O bond in 
protonated ethanol is incrementally lengthened and all other internal coordinates are 
allowed to adjust to the lowest energy configuration. The energies obtained for these 
calculations show that there is no barrier for the reverse of reaction (15), which is 
consistent with what was found using MP4 [64]. We did not find an energetically 
stable cluster between H 2 0 and ethyl cation, which is also consistent with this MP4 
study [64]. Consistent with established thermochemistry (Table 1), we find that 
reaction (15) is only 3 kcal mol -1 more endothermic than reaction (14). The rate 
constants that we obtain for reaction (15) are slightly larger than those for reaction 
(14) due to the loose transition state (large A factor) associated with reaction (15). 
The rate constants for this reaction are 12-15 orders of magnitude larger than those 
for dehydration of neat ethanol. The short lifetimes for reaction (15), r 25 o ° c = 5 x 
10 -3 s, t 5 oo °c = 2xl0 -7 s, are such that this type of reaction could also be 
important in practical pyrolysis applications. Importantly, this type of reaction 
forms a carbocation, which is a very reactive species that could be important in 
cross-linking and char formation. 

The relative rate constants that we obtain for reactions (14) and (15) are consistent 
with collisionally induced dissociation (CID) mass spectrometry experiments of 
protonated ethanol [65,66]. These experimental studies show that C 2 H^ (mlz = 29) 
is preferentially formed over H 2 0 + (m/z= 19) during collisions of protonated 
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ethanol with a collision gas. This is consistent with our calculations showing that 
reaction (15) is roughly an order of magnitude larger than reaction (14). 


3.5. Isopropanol 


Most carbohydrates have OH groups bonded to secondary carbons, and we 
investigated 1,2-dehydration of isopropanol to identify differences from ethanol, 
where the OH group is on a primary carbon. With neat isopropanol, 1,2-dehydration 
yields products analogous to those from ethanol. 


^ P H 

H2C - CH -CH 3 — 


H- 


.-pH 


hLm»..g—Q-»'J 


»«H 


’CH 3 


-v< H 

h ch 3 


+ HoO 


(16) 


The enthalpy of this reaction (A react 7/ 2 98 = 13.3 kcal mol -1 ) and the barrier to 
dehydration (AE 0 = 67.0 kcal mol -1 ) are calculated to be similar to those values for 
ethanol, Table 1. The transition state structure, Table 2, is also similar to that for 
ethanol except that the C-O bond (1.94 A) is longer and the critical dihedral angle 
(159.4°) is larger, suggesting a later transition state. The Mulliken charge distribution 
in the transition state shows a greater charge separation than in the transition state 
for ethanol, due to the ability of the secondary carbon atom to accommodate more 
positive charge. 



+0.27 

Mulliken charges in transition state 


Thus, one would expect a larger influence of added positive ions. The addition of 
Li + has a more profound effect upon the barrier for isopropanol than it does for 
ethanol, lowering the dehydration barrier 17 kcal mol -1 for isopropanol compared 
with 12 kcal mol -1 for ethanol. The result of this lowered barrier is a calculated rate 
constant for reaction (17) that is 6-8 orders of magnitude greater than the rate for 
neat isopropanol, as shown in Table 3. 


v + 

H OH 


¥ 




CH 3 


y + 

*.f H 

^CH 3 




- H//,„. 

H 1 


4= C <^^^ CH 2 CH2 + H 2 0*Li + 


(17) 


Due to the greater stability of secondary carbocations relative to primary 
carbocations, the enthalpy for the reaction (18) is roughly 12 kcal mol -1 lower 
than for ethanol. 
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^ 'r (18) 

H ^e—G ^h —- CH 3 CHCH 3 + + LiOH 
K XH 3 

However, the reaction barrier for (18) is still high, 81.5 kcal mol -1 , and the 
calculated rate constants are still 8-12 orders of magnitude smaller than those for 
reaction (17). Again, Li + greatly increases dehydration through the loss of H 2 0 
Li + , though formation of the carbocation is insignificant. 

The effect of greater stability of the carbocation is particularly impressive with 
protonated isopropanol, as is shown in Table 1 and Fig. 4. Here the formation of 
propylene, reaction (19) is about 5 kcal mol -1 more endothermic than the analogous 
channel for ethanol. 




H' 


,.-"QH 2 + 




<pn 2 + 

k=c< 


»H_ 

*ch 3 


ch 2 chch 3 + h 3 o + 


(19) 


Conversely, the energy of reaction barrier for this reaction, A£ , 0 = H-8 kcal 
mol - \ is significantly lower than the barrier for protonated ethanol, A E 0 = 20.7 kcal 
mol -1 . As was the case with ethanol, the transition state is at a lower in energy than 
the separated products and thus, we use the reaction energy as the activation energy, 
E a , instead of the energy of the transition state. The A factor was calculated with 
TST using the transition state shown in (19). The resulting rate constants are 10-14 



Fig. 4. Schematic drawing showing the energetics for the 1,2-dehydration of protonated isopropanol. The 
solid line is for reaction (19) and the dashed line is for reaction (20). The calculated relative energies, AE 0 , 
are obtained using CBS-QB3 and include zero point energy. 
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orders of magnitude greater than for neat propanol and the life times, t 2 50 °c = 3.6 s, 
T 500 °c = 4xl0 -5 s, are small enough that this reaction should be considered 
significant in pyrolysis applications. 

Significantly, greater stability of the secondary carbocation lowers the enthalpy of 
reaction (20), which is 9 kcal mol -1 more exothermic than the analogous reaction 
for protonated ethanol. 


^ pn 2 + 

™e—c^, H - 

H ^CH 3 


,oh 2 

— ch 3 chch 3 + + h 2 o 
iX ^CH 3 


( 20 ) 


This reaction (A react // 298 = 25.1 kcal mol -1 ) is now thermodynamically favored 
over reaction (19), A react // 298 = 38.1 kcal mol -1 . Additionally, the water cluster of 
the carbocation was found to be stable. Calculations of the rate constants for 
reaction (20) using the energy of reaction, 22.7 kcal mol - \ for E a and an A factor of 
10 16 s -1 yield rate constants that are 15-21 orders of magnitude greater than 1 , 2 - 
dehydrtion of neat isopropanol. The lifetimes of this reaction, r 25 o °c = 3xl0 - 7 s, 
T 500 °c = 3 x 10 -10 s, are small enough to again suggest that this type of reaction 
may be important in practical systems. Our calculations show that reaction (20) has 
a rate 5-7 orders of magnitude greater than reaction (19). This is consistent with 
CID experiments of protonated isopropanol [65]. These studies primarily show the 
formation of (mlz = 47) with no H 2 0 + (mlz= 19). In summary, secondary 

alcohols, such as are found in carbohydrates, can easily form carbocations by 1 , 2 - 
dehydration reactions at relevant temperatures. These cations may be initiators in 
char forming steps. 


3.6. Ethyleneglycol 


Ethylene glycol serves as a model system for glycol groups in carbohydrates as 
discussed above. In many respects, the 1,2-dehydration of neat ethylene glycol is 
similar to dehydration of ethanol. 







HO^ 





+ h 2 o 


( 21 ) 


As Table 1 shows, the energetics are similar, A react 7 / 298 = 7.9 kcal mol -1 , AE 0 = 
69.6 kcal mol -1 for the transition state, and as Table 2 shows, the transition states 
structures are similar. 

The similarity to ethanol disappears when ethylene glycol is complexed with Li + . 
This complex is considerably stronger than that formed with ethanol or isopropanol. 
In its complex with ethylene glycol, Li + binds to both OH groups. As a result, our 
calculated binding enthalpy for this complex is 58.4 kcal mol -1 (Table 4) while the 
binding enthalpy of Li + with ethanol and isopropanol are 38.1 and 39.8 kcal mol -1 , 
respectively. This strong bonding in the ethylene glycol results in a high activation 
energy for the loss of H 2 O Li + , reaction (22). 
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H(^ jpH 

jppG—C^ H 


V + 




H 2 Q 


HO^ 


Li + 


r 


- C = C - 


..,*H 


-Vc“! 


(22) 


In the transition state for this reaction, one of the OH groups twists away from the 
Li + and becomes planar with the ethylene group. This effectively breaks one of the 
Li-O bonds and one would expect that the binding energy of the Li + in the 
transition state to be about that of ethanol or propanol (39 kcal mol -1 ). This notion 
is supported by the similarity of the structure for this transition state (Table 2) to 
that of the transition state for lithiated ethanol. Since there is about 19 kcal mol -1 
less Li + binding energy in the transition state than in the reactant, we would expect 
the barrier for reaction (22) to be higher than the barrier for lithiated ethanol by 
about this much. The barrier for lithiated ethanol, reaction (10), is 55.3 kcal mol -1 
and we expect a barrier of about 74 kcal mol -1 for reaction (22). This is in excellent 
agreement with our calculated value, 74.2 kcal mol -1 . This high barrier results in 
rate constants for reaction (22) that are 1-2 orders of magnitude smaller than the 
dehydration rate for neat ethylene glycol. Thus, 1,2-dehydration of ethylene glycol is 
not enhanced by Li + . 

As with ethanol and isopropanol, loss of LiOH from lithiated ethylene glycol, 
reaction (23), is prohibitively high. 


A @ 

HQ'' OH © 

\ f ^OH 

(W —- H 3 G—+ Li0 H 
H H H 


(23) 


Note that for this reaction the product cation rearranges to the more stable 
protonated acetaldehyde. The enthalpy for reaction (23) is 74.5 kcal mol -1 and the 
calculated rate constants are about the same as for dehydration of neat ethylene 
glycol and thus this reaction is not likely to be important. 

Protons, unlike lithium cations facilitate dehydration of ethylene glycol (Fig. 5). 
With protonated ethylene glycol, the two OH groups do not appear to bind to the 
proton, and as a result the calculated proton affinity of ethylene glycol (187.1 kcal 
mol -1 ) is only slightly higher than the calculated proton affinity of ethanol (182.9 
kcal mol -1 ). Therefore, the calculated relative energy of the transition state is only 
27.6 kcal mol -1 above the reactant. The structure of the transition state that we 
obtain for protonated ethylene glycol is similar to that for protonated ethanol and 
protonated isopropanol as can be seen in Table 2. However, IRC calculations 
indicate that this transition state connects this transition state to protonated 
acetaldehyde and water as shown in reaction (24). 


oh® 7 oh _ 


© 


H / H 

Nc 

\ Tj 


t 


h 2 o., 

"h v r“ 

H H 


u © 


- P H 

H 3 e—CH 


© 

+ h 2 o 


(24) 
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Fig. 5. Schematic diagram of the energetics of the 1,2-dehydration of neat ethylene glycol, solid line, 
reaction (21), compared with protonated ethylene glycol, dashed line, reaction (24). The calculated relative 
energies, AE 0 , are obtained using CBS-QB3 and include zero point energy. 


The stability of the resulting oxonium cation accounts for the low enthalpy of this 
reaction, A react 7/ 2 98 = 1.6 kcal mol -1 . Since the separated products are at a lower 
energy than the transition state, we can use the TST methods outlined above to 
estimate the rate constants for reaction (24). As Table 3 shows, we obtain rate 
constants that are 13-18 orders of magnitude larger than those for neat ethylene 
glycol. Lifetimes from these rate constants, r 25 o °c=4xl0 -4 s, r 50 o °c = 6 x 10 -8 
s, are small enough that this type of reaction could be important. Oxonium ions such 
as those formed from reaction (24) could be an important source of cross-linking and 
char formation in carbohydrate pyrolysis conducted in the presence of acid. 


3. 7. Implications for carbohydrate pyrolysis 

Simple 1,2-dehydrations of all of the neat alcohols studied consistently have 
barriers of approximately 67-69 kcal mol -1 . As Table 1 shows this barrier is 
independent of size of the alcohol and is the same for glycols groups. Thus, it is 
reasonable to assume that this barrier is the same for the carbohydrates found in 
biomass. We have computed the barrier for dehydration in a simple sugar, 
levoglucosan, to test this hypothesis. 
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+ H 2 0 


( 25 ) 


Levoglucosan 


1 


Since this molecule is a bicyclic species there are a limited number of conformers, 
which makes the computation easier. It is impractical to conduct CBS-QB3 
calculations of this molecule, but the computed barrier for dehydration using 
B3LYP is possible on a large computer. Our calculations show that the B3LYP 
barrier, 68 kcal mol -1 , is consistent with the calculated barriers for other alcohols, 
as is shown in Table 1 and it, therefore, seems reasonable to assume that the actual 
barrier of 1,2-dehydration for levoglucosan and other sugars is 67-69 kcal mol -1 . 


4. Conclusions 

From our investigation of the 1,2-dehydration of alcohols, we can draw several 
conclusions concerning the importance of these reactions, (i) In general, the addition 
of alkali metal cations and protons dramatically increases the rates of these 
reactions. With protonation, the calculated rates are high enough that this process 
appears likely on a time-scale of typical pyrolysis or torrefaction processes. The 
increase in rates is primarily due to a decrease in the potential barrier, which is 
caused by stabilization of transition state. The trend for this effect in terms of 
decreasing energy of the reaction barrier is neat alcohol <Na + <Li + <H + . The 
trend appears to be related to the size of the cation and one would expect that of the 
two important alkali metals in biomass, K + would have a smaller impact than Na + . 
(ii) This decrease in the potential barrier becomes larger for a more highly 
substituted alcohol. This is likely due to greater accommodation of charge by 
secondary carbons relative to primary carbons, (iii) Greater substitution of the 
alcohol also leads to the formation of carbocations and oxonium as products. These 
highly reactive products may serve as cross-linking agents, which can lead to the 
formation of char during pyrolysis, (iv) Strong metal/glycol bonding slows 1,2- 
dehydration, because the glycol is stabilized relative to the transition state. 

Consider the pyrolysis protonated cellulose. The following shows a potential 
reaction of protonated cellulose under pyrolysis or torrefaction conditions. 



(26) 



If our studies of model compounds are applicable, it is reasonable to expect that 
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the first step is fascile under typical conditions. However, this picture may be 
complicated by the hydrogen bonding found in solid cellulose. The effect of 
additional hydrogen bonding can and should be investigated in the future. In any 
event, the reactivity of the resulting oxonium ion makes the subsequent reaction 
likely. 
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